ELSEVIER 


Available  online  at  www.sciencedirect.com 


SCIENCE 


DIRECT® 


Journal  of  Power  Sources  156  (2006)  92-99 


www.elsevier.com /locate /jpowsour 


Study  of  gas  pressure  and  flow  rate  influences  on  a  500  W  PEM  fuel  cell, 

thanks  to  the  experimental  design  methodology 

B.  Wahdamea>*,  D.  Candussob,  Jean-Marie  Kauffmanna 

a  Laboratory  of  Electrical  Engineering  and  Systems  (L2ES)  -  EA  3898,  University  of  Technology  Belfort-Montbeliard  (UTBM)  &  University  of 

Franche-Comte  (UFC),  L2ES-UTBM  Bdtiment  F,  Technopole,  Rue  Thierry  Mieg,  F  90010  Belfort,  France 
b  The  French  National  Institute  for  Transport  and  Safety  Research  (INRETS),  Laboratory  of  New  Technologies  (LTN), 

2  Avenue  du  General  Malleret-Joinville,  F  94114  Arcueil  Cedex,  France 

Available  online  25  October  2005 


Abstract 

The  behaviour  of  a  500  W  PEM  fuel  cell  stack,  fed  by  pure  hydrogen  and  humidified  compressed  air,  is  currently  investigated  on  the  fuel  cell 
test  platform  of  Belfort. 

In  this  paper,  the  influences  on  fuel  cell  performance  of  gas  pressure  and  flow  rate  parameters  are  studied.  The  fuel  cell  is  operated  in  the  pressure 
regulation  mode:  the  gas  flow  rates  are  regulated  thanks  to  mass  flow  controllers  placed  upstream  of  the  stack  and  the  gas  pressures  at  stack  inlets 
are  controlled  by  regulation  valves  located  downstream  of  the  stack.  The  choice  of  the  various  tests  to  perform  is  made  thanks  to  experimental 
design  methodology,  which  is  a  suitable  technique  to  characterise,  analyse  and  to  improve  a  complex  system  such  as  a  fuel  cell  generator.  In  this 
study,  the  four  physical  factors  considered  are  both  hydrogen/air  pressures  and  anode/cathode  flow  rates.  Each  factor  has  two  levels,  leading  to 
a  full  factorial  design  requiring  16  experiments  (16  current- voltage  curves).  The  test  bench  developed  at  the  laboratory  allows  setting  the  other 
factors  (for  instance:  stack  temperature,  relative  humidity  and  dew  point  temperature  of  the  air  at  stack  inlet)  at  fixed  values.  The  test  responses  are 
the  maximal  output  power  and  the  efficiency  computed  for  this  power.  Statistical  sensitivity  analyses  (ANOVA  analyses)  are  used  to  compute  the 
effects  and  the  contributions  of  the  various  factors  to  the  fuel  cell  maximal  power.  The  use  of  fractional  designs  shows  also  how  it  is  possible  to 
reduce  the  number  of  experiments.  Some  graphic  representations  are  employed  in  order  to  display  the  results  of  the  statistical  analyses  made  for 
different  current  values. 

©  2005  Elsevier  B.Y.  All  rights  reserved. 
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1.  Introduction 

Proton  exchange  (or  polymer  electrolyte)  membrane  fuel  cell 
(PEMFC)  is  considered  as  being  one  of  the  most  promising  tech¬ 
nologies  able  to  produce  efficient  and  environmentally  friendly 
energy  for  various  applications.  More  and  more  attention  is  paid 
to  PEMFC  for  powering  electric  vehicles  especially  because 
of  its  low  temperature  operation  and  its  high  power  density. 
Nevertheless,  PEMFC  is  not  yet  an  economical  solution  for  the 
transportation  market  area  and  some  research  have  still  to  be 
led  in  order  to  reduce  the  FC  generators  costs,  to  improve  their 
reliability  as  well  as  to  find  some  solutions  for  the  hydrogen 
distribution. 
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The  F2ES  laboratory  based  on  the  fuel  cell  test  platform  of 
Belfort  is  currently  carrying  research  on  a  500  W  PEMFC  stack. 
The  aim  is  to  ensure  for  the  FC  proper  and  optimal  operating 
conditions,  leading  to  high  power  or  efficiency  delivery.  This 
task  is  not  easy  especially  because  of  the  high  number  of  control 
parameters  such  as:  FC  temperature,  gas  pressure,  gas  flow,  rel¬ 
ative  humidity  of  hydrogen  and  air,  profile  of  the  load  current, 
etc.  All  these  parameters  have  strong  impacts  on  the  FC  voltage 
and  are  related  among  themselves  by  nonlinear  relations,  diffi¬ 
cult  to  be  modelled.  Their  values  depend  on  the  technological 
choices  made  for  the  ancillaries  (for  instance,  the  air  pressure 
and  flow  levels  are  linked  to  the  selected  type  of  compressor  and 
to  its  characteristics). 

A  large  number  of  experimental  tests  are  often  needed  to 
correctly  analyse  the  performances  of  a  given  FC  system  or 
to  identify  the  parameters  of  a  physical  model.  The  design  of 
experiment  (DOE)  method  can  be  used  in  order  to  evaluate  the 
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Table  1 


Technical  specifications  of  the  fuel  cell 


Number  of  cells 

20 

Cell  area  (cm2) 

100 

Operating  temperature  (°C) 

20-65 

Operating  pressure 

Maximum  1.5barabs  (0.5  barreiatiVe) 

Nominal  power  output  (W) 

500 

Operating  power  (W) 

0-800 

Media  inlet 

Cathode 

Humidified  air 

Anode 

Pure  hydrogen  (dry) 

Cooling 

Demineralised  water 

respective  impacts  of  the  physical  control  parameters  on  the  FC 
operation.  The  DOE  method  dates  in  fact  back  to  the  beginning 
of  the  last  century  with  the  work  from  Fisher  (1925).  The  first 
users  of  these  methods  were  agronomists  who  quickly  included 
the  interest  of  the  experimental  designs,  in  particular  the  pos¬ 
sible  reduction  of  the  number  of  tests  when  many  parameters 
are  studied.  In  the  1960s,  the  DOE  brought  many  innovations, 
mainly  thanks  to  the  work  done  by  Dr.  Genichi  Taguchi  [1-3]. 
This  methodology  is  sometimes  used  by  chemists  [6]  in  order  to 
characterise  the  fuel  cell  and  its  materials,  to  determine  the  most 
significant  parameters,  or  to  highlight  the  possible  correlations 
between  these  ones. 

In  our  work,  DOE  tests  have  been  conducted  with  two  objec¬ 
tives:  to  obtain  the  maximum  of  the  power  produced  by  the  fuel 
cell  and  to  reach  the  maximum  of  the  fuel  cell  efficiency  for 
given  ranges  of  air/hydrogen  flows  and  pressures.  Experiment 
results  are  given  and  analysed  in  this  paper. 

2.  Experimental  setup 

The  PEMFC  used  in  this  study  has  been  assembled  with  a 
sprayed  catalyst  layer  (Gore  MESGA  Primea  Series  5510)  and 
graphite  distribution  plates.  Some  data  on  the  fuel  cell  are  given 
in  Table  1 . 

In  order  to  achieve  the  objectives  of  this  study,  we  have 
adapted  the  DOE  method  to  the  available  testbench  in  our  labora¬ 
tory  (Fig.  1)  [7,8].  The  500  W  stack  was  operated  at  atmospheric 
pressure  with  humidified  air  (dew  point  temperature  of  25  °C 


Table  2 


The  levels  used  for  the  various  experiments 


Factors 

Levels 

Minimum  (/=  —  1) 

Maximum  (i  =  +1) 

A:  Ph2 

—  1:  1.3  bar 

1:  1.4  bar 

B-  Pair 

—  1:  1.3  bar 

1:  1.4  bar 

C :  Dh2 

—  1:  5.5  nlmin-1 

1:  8.5  nlmin-1 

D  '-  Dair 

—  1:  26.3  nlmin-1 

1:  39.6 nlmin-1 

and  relative  humidity  of  75%  at  stack  inlet)  and  dry  hydrogen. 
The  air  and  hydrogen  flows  were  controlled  by  flow  regulators 
placed  upstream  of  the  stack.  The  pressures  at  stack  inlets  were 
controlled  thanks  to  back  pressure  valves  located  downstream 
of  the  stack.  The  stack  temperature  was  regulated  at  55  °C. 

During  each  one  of  the  16  tests,  hydrogen  and  air  flow  rates 
were  kept  constant.  For  example,  during  the  experiment  pre¬ 
sented  below  (Fig.  2),  their  values  were,  respectively,  8.5  and 
39.6  nl  min-1.  The  hydrogen  and  air  pressures  were  controlled 
progressively  as  the  current  was  incremented  by  5  or  2  A  steps. 
The  output  power  increased  progressively  and  so  did  the  volt¬ 
age  drop  of  the  stack.  When  any  one  of  the  20  cell  voltages  was 
reaching  the  threshold  of  380  mV,  the  test  was  stopped. 

3.  Full  factorial  design 

A  design  of  experiment  (DOE)  is  a  structured,  organised 
method  to  determine  the  relationship  between  the  factors  (X) 
affecting  a  process  and  the  output  of  that  process  (I7)-  Here,  the 
experimental  design  factors  considered  are:  A:  hydrogen  pres¬ 
sure  (Eh?)>  B\  air  pressure  (Fair),  C :  hydrogen  flow  rate  (Dh?), 
and  D  which  is  the  air  flow  rate  (Dair).  The  response  Y  is  either 
the  maximum  FC  output  power  or  the  efficiency  computed  for 
this  power.  The  levels  adopted  for  the  factors  are  summed  up  in 
Table  2. 

The  levels  were  chosen  taking  into  account  the  limits  of  the 
FC  and  also  the  ones  of  the  testbench.  Especially  the  limited 
pressure  operating  range  specified  by  the  FC  manufacturer  was 
restrictive  for  our  tests.  On  the  other  hand,  it  could  be  considered 
that  the  small  amplitude  of  pressure  and  gas  flow  parameters  let 
the  FC  operate  in  a  homogeneous  physical  domain.  Therefore 
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Fig.  1.  Simplified  scheme  of  the  fuel  cell  testbench. 
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Fig.  2.  (a)  Evolution  of  the  FC  current,  (b)  evolution  of  the  20  cell  voltages,  (c)  evolution  of  the  FC  stack  power  and  (d)  evolution  of  the  gas  pressures  at  stack  inlets. 


the  hypothesis  of  linearity  for  these  factors  could  be  considered. 
This  hypothesis  had  to  be  verified  afterwards. 

DOE  can  be  very  simple  or  very  complicated  according  to 
the  number  of  studied  factors  and  formulated  assumptions.  The 
full  factorial  design  of  four  presumed  linear  factors  includes  24 
experiments  (four  factors  with  two  levels  each)  (Table  3)  [1,2,4]. 


Table  3 

Full  factorial  format  for  four  linear  factors 


Test  no. 

A 

B 

C 

D 

ABCD 

1 

2 

+ 

— 

— 

— 

+ 

3 

— 

+ 

— 

— 

— 

4 

+ 

+ 

— 

— 

+ 

5 

— 

— 

+ 

— 

— 

6 

+ 

— 

+ 

— 

+ 

7 

— 

+ 

+ 

— 

+ 

8 

+ 

+ 

+ 

— 

— 

9 

— 

— 

— 

+ 

— 

10 

+ 

— 

— 

+ 

+ 

11 

— 

+ 

— 

+ 

+ 

12 

+ 

+ 

— 

+ 

— 

13 

— 

— 

+ 

+ 

+ 

14 

+ 

— 

+ 

+ 

— 

15 

— 

+ 

+ 

+ 

— 

16 

+ 

+ 

+ 

+ 

+ 

For  the  effect  calculation,  we  have  used  the  method  of  Yates. 
The  notation  used  in  Yates  method  is  very  practical  if  the  studied 
factors  have  two  levels.  For  instance,  the  effect  A  on  level  i  is 
computed  as  follows: 

Em  =  hi  ~  Y  (1) 

/V 

where  Yai  is  the  average  of  the  responses  when  the  factor  A 

/V 

is  on  level  i  and  Y  is  the  average  of  the  responses  for  all  the 
experiments. 

In  order  to  compare  the  eight  effects  of  the  four  factors,  we 
represent  these  effects  on  a  graph  (Fig.  3)  that  is  called:  graph 
of  the  average  effects.  The  maximal  powers  reached  for  the  16 
experiments  are  plotted  on  the  four  graphs  below,  as  well  as  the 
general  average  of  these  167  values  and  the  average  of  the  stack 

/V 

powers  when  the  considered  X  factor  is  on  levels  + 1/  —  l(Yxi). 
On  the  four  graphs  below,  the  greater  the  slope  of  the  segment 

/V  V\ 

[y^(_  i)  the  larger  is  the  influence  of  the  X  factor  [1,2]. 

The  analysis  of  variance  (ANOVA)  is  a  commonly  used  tool 
to  study  and  to  estimate  the  factor  influences  over  a  process.  A 
first  ANOVA  analysis  is  performed  for  the  16  reached  maximal 
powers.  Only  first  degree  interactions  (AB  for  instance)  are  taken 
into  account.  The  analysis  results  are  summed  up  in  Table  4. 

The  ANOVA  table  has  seven  columns.  The  first  one  shows 
the  source  of  the  variability.  The  second  one  shows  the  sum  of 
squares  due  to  each  source.  The  third  one  includes  the  percent- 
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Influence  of  A  (H2  pressure)  on  Y 


(a)  Variable  A  (H2  pressure)  (-1  =1.3  bars)  (1  =1.4  bars) 


Influence  of  B  (air  pressure)  on  Y 
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(c)  Variable  C  (H2  flow  rate)  (-1  =  5.5  nl/min)  (1  =  8.3  nl/min) 


Influence  of  D  (air  flow  rate)  on  Y 


(d)  Variable  D  (air  flow  rate)  (-1  =  26.3  nl/min)  (1  =  39.6  nl/min) 


Fig.  3.  Graphs  of  the  mean  effects  for  the  four  factors. 


age  of  contribution  for  each  factor.  The  fourth  one  represents 
the  degrees  of  freedom  (d.f.)  associated  with  each  source.  The 
fifth  one  contains  the  mean  squares,  which  is  the  ratio:  sum  of 
squares/d.f.  The  sixth  one  shows  the  Fisher  statistics  (fs),  which 
is  the  ratio  of  the  mean  squares.  The  seventh  one  shows  the 
p- value  for  the  Fisher  statistics.  The  choice  of  a  limit  for  the 
jp-value,  in  order  to  determine  whether  or  not  a  result  is  “sta¬ 
tistically  significant”,  is  left  to  the  researcher.  It  is  common  to 


declare  a  result  significant  if  the  p-  value  is  less  than  0.05  or  0.01 

[1,2]. 

The  ANOVA  of  Table  4  shows  that  the  air  flow  rate  (factor  D) 
is  the  most  important  factor  with  a  contribution  equal  to  88.5% 
of  the  total  variance.  The  effect  of  the  factor  C  is  the  second 
most  significant. 

The  reversible  cell  efficiency  r/rzv  is  the  ratio  between  the 
electrical  power  output  (Pel)  and  the  product  of  the  total  enthalpy 


Table  4 

ANOVA  for  the  design  of  experiments  LI 6  with  maximal  powers 


Source 

Sum  of  squares 

Percentage 

d.f. 

Mean  square 

Fisher  statistics 

P 

A  ( Pu2 ) 

749.4 

0.8 

1 

749.4 

2.99 

0.144 

B  (Pair) 

199.5 

0.21 

1 

199.5 

0.8 

0.413 

C(Dh2) 

5531.6 

6 

1 

5531.6 

22.1 

0.005 

D  (f^air) 

82010 

88.5 

1 

82010 

327 

0 

AB 

21.4 

2E— 3 

1 

21.4 

0.09 

0.781 

AC 

78.8 

5E-3 

1 

78.8 

0.31 

0.599 

AD 

8.3 

0 

1 

8.3 

0.03 

0.863 

BC 

28.9 

0.03 

1 

28.9 

0.12 

0.747 

BD 

213.9 

0.23 

1 

213.9 

0.85 

0.397 

CD 

2512.5 

2.71 

1 

2512.5 

10 

0.024 

Error 

1252.1 

1.35 

5 

250.4 

Total 

92607 

15 
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Table  5 


ANOVA  for  DOE  L16  with  efficiencies  at  maximal  powers 


Source 

Sum  of  squares 

Percentage 

d.f. 

Mean  square 

Fisher  statistics 

P 

A  (Ph2) 

6.13 

0.6 

1 

6.126 

3.02 

0.142 

B  (Pair) 

1.05 

0.1 

1 

1.051 

0.52 

0.503 

C(Dh2) 

487.31 

47 

1 

487.306 

240.3 

0 

D  (Da ir) 

527.85 

50.9 

1 

527.851 

260.3 

0 

AB 

0.05 

0 

1 

0.051 

0.02 

0.880 

AC 

1.5 

0.14 

1 

1.501 

0.74 

0.428 

AD 

0.46 

0.04 

1 

0.456 

0.22 

0.655 

BC 

0.02 

0 

1 

0.016 

0.01 

0.933 

BD 

1.76 

0.17 

1 

1.756 

0.87 

0.394 

CD 

0.86 

0.08 

1 

0.856 

0.42 

0.544 

Error 

10.14 

0.97 

5 

2.028 

Total 

1037.1 

15 

Table  6 

ANOVA  table  for  the  fractional  factorial  design  (FC  maximal  power  is  the  output) 

Source 

Sum  of  squares 

Percentage 

d.f. 

Mean  square 

Fisher  statistics 

P 

A  (Ph9) 

247.5 

0.52 

1 

247.5 

0.48 

0.538 

B  (Pair) 

52.5 

0.11 

1 

52.5 

0.1 

0.770 

C(P>h2) 

2397.8 

5 

1 

2397.8 

4.64 

0.120 

D  (D air) 

42997.8 

91 

1 

42997.8 

83.26 

0.003 

of  the  reaction  ( AHf )  with  the  hydrogen  molar  flow  rate  (Vh2), 

generator  D  =  ±ABC  and  ABCD  =  I  (identity)  (Table  3).  In  this 

as  expressed  by  formula  (2)  [5]. 

way,  a  fractional  factorial  design  is  generated  [2]. 

P  i 

ANOVA  is  carried  out  for  the  eight  powers;  the  interactions 

7?  rev  — 

r  el 

(2) 

are  not  taken  into  account  (Table  6). 

-Vh2  x  A  iff 

The  results  obtained  by 

means  of  the  fractional  design  are 

The  efficiency  is  computed  for  each  one  of  the  16  maximal 

close  to  the  ones  of  the  full  factorial  design;  a  very  significant 

electrical  powers  and  a  new  ANOVA  analysis  is  performed. 

effect  of  the  parameter  D  (air  flow  Da ir)  is  found;  it  represents 

In  Table  5,  we  can  see  that  the  efficiency  of  the  500  W  FC 

91  %  of  the  total  of  contributions.  In  the  same  way,  the  efficiency 

shows  an  important  dependence  on  both  factor  D ,  with  50.9% 

study  with  a  fractional  factorial  design  reveals,  as  expected  that 

of  the  total  variance,  and  factor  C  with  47%. 

hydrogen  and  air  flow  rate  influences  are  dominant  over  the  fuel 

cell  stack  efficiency  (Table  7). 

4.  Fractional  factorial  design 

The  1 6  experiments  and  the  ANOVA  analyses  performed  with 

a  full  factorial  design  show  that  an  increase  of  air  flow  is  the  best 

A  fractional  factorial  DOE  includes  selected  combinations  of 

way  to  obtain  the  highest  power  and  efficiency,  for  the  range  of 

factors  and  levels.  It  is  a  carefully  prescribed  and  representative 

flows  and  pressures  considered.  The  use  of  the  fractional  design 

subset  of  a  full  factorial  design.  A  fractional  factorial  DOE  is 

leads  to  the  same  conclusion,  but  only  eight  instead  of  16  exper- 

useful  when  the  number  of  potential  factors  is  relatively  large 

iments  are  needed. 

because  it  reduces  the  total  number  of  experiments  required.  A 

fractional  design  has  to  check  the  following  properties:  the  first 

5.  Experimental  result  analysis  and  discussion 

one  is  orthogonality  and  the  second  one  is  linked  to  the  degrees 

of  freedom. 

5.7.  The  effects  of  pressure 

The  full  factorial  design  required  2k  experiments  where  k  is 

the  number  of  two  level  factors.  From  a  full  factorial  design,  it 

It  might  be  wondering  that  the  effects  of  pressures  are  so 

is  possible  to  reduce  the  runs  to  8  (=. 

24  thanks  to  the  design 

low,  even  if  the  pressure  range  considered  is  small  (0.1  bar)  as  it 

Table  7 

ANOVA  table  for  the  fractional  factorial  design  (FC  efficiency  computed  for  the  maximal  power  is  the  output) 

Source 

Sum  of  squares 

Percentage 

d.f. 

Mean  square 

Fisher  statistics 

P 

A  (Ph2) 

1.62 

0.3 

1 

1.62 

8.53 

0.061 

B  (Pair) 

0.605 

0.1 

1 

0.605 

3.18 

0.172 

C(P>h2) 

262.205 

48.7 

1 

262.205 

1380 

0 

D  (Pair) 

273.78 

50.8 

1 

273.78 

1440.9 

0 
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Table  8 


ANOVA  for  DOE  L16  (FC  power  computed  for  a  20  A  current  is  the  output) 


Source 

Sum  of  squares 

Percentage 

d.f. 

Mean  square 

Fisher  statistics 

P 

A  (Ph2  ) 

30.582 

8.95 

1 

30.582 

2.14 

0.203 

B  (P air) 

51.407 

15.05 

1 

51.407 

3.6 

0.116 

C(Dh2) 

3.535 

1.035 

1 

3.535 

0.25 

0.639 

19  iP air) 

165.897 

48.58 

1 

165.897 

11.63 

0.019 

AB 

5.429 

1.597 

1 

5.429 

0.38 

0.564 

AC 

5.382 

1.576 

1 

5.382 

0.38 

0.565 

AD 

2.823 

0.826 

1 

2.823 

0.2 

0.675 

BC 

0.961 

0.28 

1 

0.961 

0.07 

0.805 

BD 

0 

0 

1 

0 

0 

0.996 

CD 

4.121 

1.206 

1 

4.121 

0.29 

0.613 

Error 

71.305 

20.88 

5 

14.261 

Total 

341.443 

15 

has  already  been  noticed.  Indeed,  the  Nernst  equation  and  also 
the  expression  of  the  activation  overvoltages  clearly  show  that 
increasing  the  pressures  leads  to  a  higher  fuel  cell  voltage.  The 
Nernst  equation  can  be  written  as  follows  [5]: 


E  =  E°  + 


RT 

2F 


In 


(  Pu2  x 

\  PH20  J 


where  E°  (about  1.2  V)  is  the  limit  voltage  at  standard  pressure, 
Po2  the  partial  pressure  of  oxygen  (bar),  Pu2o  the  partial  pres¬ 
sure  of  water  (bar)  and  Pn2  the  partial  pressure  of  hydrogen 
(bar). 

The  expression  for  the  activation  voltages  can  be  derived  from 
the  Butler- Volmer  equation: 

rizct  =  +  frT  +  foT In  +  p4T  In (i) 


+  05  In 


where  i  is  the  current  density  (Am  2),  Po  the  reference  pressure 
( Po  =  1  atm). 

Each  one  of  the  16  tests  performed  can  be  used  to  record  a 
polarisation  curve.  In  Fig.  4,  two  polarisation  curves  are  plotted 
in  the  case  of  two  different  pressure  levels. 


Fig.  4.  Polarisation  curves  recorded  for  two  hydrogen  and  air  pressures. 


The  gas  pressure  increase  from  1.3  to  1.4  bar  leads  to  an 
improvement  of  the  FC  performances.  This  contribution  of  the 
pressures  to  the  FC  voltage  output  (and  thus  to  the  power  output) 
was  not  noticed  in  the  results  of  the  first  experimental  design, 
where  we  were  only  interested  in  the  maximum  power  reached. 
A  new  study  has  been  carried  out  for  another  operation  point  on 
the  polarisation  curve  (for  the  same  current  value  of  20  A  instead 
of  the  different  current  values  corresponding  to  the  maximal 
powers  reached)  with  the  aim  to  highlight  the  incidence  of  the 
pressures  on  the  power  delivered  by  the  fuel  cell.  The  results 
obtained  by  this  procedure  are  shown  in  Table  8. 

The  effects  of  the  factors  D  (Da ir),  B  (Pa ir)  and  A  (Ph2)  are 
more  significant  but  they  do  not  explain  the  total  contribution.  In 
fact,  there  is  a  large  error  due  to  other  influential  parameters  not 
taken  into  account  [2] .  The  analysis  of  the  variance  performed 
for  the  efficiency  of  the  fuel  cell  shows  the  very  important  effect 
of  the  hydrogen  flow  (Table  9). 

5.2.  The  effects  of  flow  rate 

It  is  well-known  that  water  management  is  of  great  impor¬ 
tance  for  the  PEMFC  correct  operation.  Not  only  the  membrane 
needs  a  certain  moisture  to  reach  a  good  conductivity;  the  FC 
performances  depend  also  on  the  amount  of  water  in  the  backing 
layers.  In  the  case  of  a  flooding,  the  catalytic  sites  in-depth  are 
not  fed  by  oxygen  (“too  wet”  case).  When  the  pores  become  too 
dry,  the  proton  supply  is  not  ensured  (“too  dry”  case).  The  two 
cases  lead  to  a  reduction  in  the  number  of  active  sites  and  thus 
to  a  performance  decrease.  For  the  correct  operation  of  our  FC, 
the  air  humidity  must  be  controlled.  The  studied  FC  needs  quite 
high  stoichiometry  rates  to  operate  properly.  It  seems  that  the  air 
flow  contributes  actively  to  the  evacuation  of  the  water  produced 
in  the  FC.  The  flow  rates  of  gases  have  a  strong  influence  on  the 
moisture  in  our  FC  and  especially  on  the  draining  away  of  the 
water  from  plate  channels  and  gas  diffusion  layers. 

5.3.  Additional  results 

The  tools  developed  for  one  particular  current  value  can  be 
used  in  order  to  understand  the  contribution  and  the  influence  of 
each  parameter  over  the  FC  voltage,  along  the  curve  of  polar- 
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Table  9 


ANOVA  for  DOE  L16  (FC  voltage  efficiency  computed  for  a  20  A  current  is  the  output) 


Source 

Sum  of  squares 

Percentage 

d.f. 

Mean  square 

Fisher  statistics 

P 

A  (Pi h) 

0.176 

0.039 

1 

0.176 

1.51 

0.274 

B  (Pair) 

0.336 

0.075 

1 

0.336 

2.87 

0.150 

C(Dh2) 

444.577 

99.42 

1 

444.577 

3795.5 

0 

D  (Dak) 

1.221 

0.273 

1 

1.221 

10.42 

0.023 

AB 

0.062 

0.014 

1 

0.062 

0.53 

0.497 

AC 

0.009 

0.002 

1 

0.009 

0.08 

0.792 

AD 

0.027 

0.006 

1 

0.027 

0.23 

0.650 

BC 

0.002 

0 

1 

0.002 

0.02 

0.900 

BD 

0.006 

0 

1 

0.006 

0.05 

0.835 

CD 

0.16 

0.036 

1 

0.16 

1.37 

0.295 

Error 

0.586 

0.001 

5 

0.117 

Total 

447.163 

15 

Current  [A] 

Fig.  5.  Contribution  of  the  four  studied  parameters  to  the  FC  voltage. 

isation.  For  instance,  ANOVA  analyses  are  performed  for  six 
various  FC  currents  in  the  case  of  the  FC  voltage  study.  The 
contributions  are  plotted  on  the  two  graphs  of  Fig.  5  (interac¬ 
tions  are  not  taken  into  account). 

6.  Conclusion 

The  design  of  experiment  method  has  been  used  in  order 
to  compute  the  effects  of  four  factors  (hydrogen  pressure,  air 
pressure,  hydrogen  flow  rate  and  air  flow  rate)  over  the  max¬ 
imum  power  delivered  by  the  FC  and  also  over  its  efficiency 
at  this  operation  conditions.  First,  two  studies  were  executed:  a 
full  factorial  design  (16  experiments)  and  a  fractional  one  (eight 
experiments).  In  both  cases,  we  arrived  to  the  same  conclusion 
concerning  the  prevalent  effects  of  the  maximal  air  flow  rate 
and  hydrogen  flow  rate,  respectively,  on  the  maximal  power 
and  efficiency  obtained.  The  studied  FC  needs  relative  high 
stoichiometry  rates  to  deliver  its  maximal  power.  The  major 
importance  of  the  air  flow  highlighted  in  the  study  could  explain 


the  fact  that  FC  manufacturers  choose  to  develop  “air  breathing” 
technologies  (or  “convection  style  stacks”).  Their  operation  type 
is  very  simple  as  the  opened  cathode  can  be  fed  by  a  blower.  As 
there  is  no  need  for  a  compressor,  the  ancillary  consumption 
is  low.  This  kind  of  FC  seems  well  adapted  for  low  powers 
(<1  kW).  In  the  second  part  of  the  study,  it  was  noticed  that  the 
pressures  effects  were  very  low,  which  seemed  to  be  in  contra¬ 
diction  with  some  simple  physical  considerations.  Therefore  a 
new  study  was  done:  we  were  no  longer  interested  in  the  maxi¬ 
mal  power  reached  but  in  the  power  supplied  at  a  given  current. 
The  results  were  appreciably  different  from  the  initial  work  inso¬ 
far  as  it  became  possible  to  highlight  the  influences  of  air  and 
hydrogen  pressures  on  the  power  delivered. 

Future  works  will  concern  the  study  of  other  physical  factors 
(for  instance:  relative  humidity  and  dew  point  temperature  of 
gases  at  stack  inlets).  The  methodologies  initially  developed  for 
the  500  W  stack  will  be  used  in  the  framework  of  other  projects 
concerning  larger  PEM  stacks  (5  kW  for  instance).  It  could  be 
also  interesting  to  consider  a  physical  model  and  to  use  the  results 
of  the  experiment  analyses  with  the  aim  at  evaluating  the  model 
and  its  degree  of  complexity  for  a  given  application. 
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